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Abstract 

The theory of resonant generation of nonground-state Bose-Einstein condensates 
is extended to Bose-condensed systems at finite temperature. The generahzation is 
based on the notion of representative statistical ensembles for Bose systems with broken 
global gauge symmetry. Self-consistent equations are derived describing an arbitrary 
nonequilibrium nonuniform Bose system. The notion of finite-temperature topological 
coherent modes, coexisting with a cloud of noncondensed atoms, is introduced. It is 
shown that resonant generation of these modes is feasible for a gas of trapped Bose 
atoms at finite temperature. 
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1 Introduction 

Statistical systems with Bose-Einstein condensate exhibit a variety of very interesting phe- 
nomena, both equilibrium as well as nonequilibrium (see, e.g., book [1] and review articles 
[2-5]). One such a nontrivial effect is the possibility of creating nonground-state Bose- 
Einstein condensates of trapped atoms, as was advanced in Ref. [6]. The properties of these 
nonground-state condensates have been studied in a number of papers, for instance in Refs. 
[6-19], and are recently summarized in survey [20]. All these works [6-20] deal solely with 
the case of trapped atoms forming a dilute gas at zero-temperature, when all atoms can be 
condensed, so that there is no noncondensed atoms. However the latter are always present 
in real experiments at finite temperature. Interactions between atoms also produce an ad- 
mixture of noncondensed particles. Then the question arises whether the nonground-state 
condensates can be generated, when the trapped atoms do not form a purely condensed 
system, but there is also a cloud of noncondensed atoms. 
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To answer the above question, it is necessary to be based on a rehable general theory 
of nonequilibrium nonuniform Bose-condensed systems. The description of such systems is 
commonly done by employing the Bogolubov ideas [21-24], when the global gauge symmetry 
is broken by means of the Bogolubov shift. However, Hohenberg and Martin [25], first, 
emphasized that all theories of Bose systems with broken gauge symmetry suffer from an 
internal inconsistency, being either nonconserving or displaying an unphysical gap in the 
collective spectrum. In both these cases, such theories do not correspond to stable statistical 
systems. A detailed discussion of this problem has recently been done by Andersen [3]. 

To overcome the standard deficiency of theories with broken gauge symmetry, it is neces- 
sary to invoke the notion of representative ensembles [26] . This notion is strictly formulated 
in Sec. 2 for arbitrary statistical systems. The appropriate representative ensemble for a 
Bose system with broken gauge symmetry is constructed in Sec. 3. The resulting evolu- 
tion equations arc derived from the general variational principle for the extremum of action 
functional, which makes all equations self-consistent and the theory conserving and gapless, 
as is shown in Sec. 4. The derived evolution equation for the condensate wave function is 
analysed in Sec. 5, where topological coherent modes at finite temperature are defined. Res- 
onant generation of these coherent modes, corresponding to nonground- state condensates, is 
demonstrated to be feasible even in the presence of a substantial admixture of noncondensed 
atoms. 

Throughout the paper, the system of units is used, where the Planck and Boltzmann 
constants are set to unity, h — 1, ks — 1- 

2 Representative statistical ensembles 

Statistical systems are characterized by statistical ensembles, which means the following. 
First, one has to specify a space J-' of microstates spanning all admissible quantum or dy- 
namical states of the given system. Second, a statistical operator p(0) at the initial time t — 
has to be fixed. Third, temporal evolution of either the statistical operator or of observable 
quantities is to be defined, which is given by defining the action of the time derivative d/dt. 
Thus, a general nonequilibrium statistical ensemble can be denoted as a triplet 

{^.P(0,.|) . (1) 

For stationary or, in particular, equilibrium statistical systems, when their evolution is 
either trivial or absent, an equilibrium statistical ensemble is a pair {J^,p}, where p{t) — 
P(0) = p. 

The temporal evolution equation of operators in the Heisenberg representation can be 
symbolized with the help of the evolution operator U (t) , so that the evolution of an operator 
C{t) is given by the relation 

C{t) = U+{t)C{0)U{t) . 
Then the statistical average of this operator is 

< C{t) > = Tr p{0)C{t) = Tr p{t)C{0) , 

where the trace operator is accomplished over the given state of microstates J-'. 
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Each statistical system is characterized by a set of dynamical variables. The latter, 
keeping in mind quantum statistical systems, can be called field variables, whose example 
is a set ip{x,t) = [iljj{x,t)] of field operators iljj{x,t), with an index j enumerating the set 
members. Here t is time and x is a collection of all other variables, which, e.g., could be 
spatial coordinates or momenta. 

The system energy operator is given by a Hamiltonian H[ijj], which is a functional of the 
field variables -0- The related Lagrangian is 

L[i;] = J i;\x,t) t ^ ^{x,t) dx - H[i;] , (2) 

where the set ip — [ipj] can be treated as a column. 

Strictly speaking, for correctly defining a statistical system, it is not always sufficient to 
fix just a Hamiltonian H[ip] or a Lagrangian L[iIj], but it is necessary to formulate additional 
conditions or constraints for making the system uniquely defined. Suppose there is a family 
{C'j['V']} of self-adjoint operators Ci[ip] = C^[ip], which will be called condition operators, 
when the required additional constraints, imposed on the system, are formulated as statistical 
conditions 

Q = < > , (3) 

defined as the statistical averages of the condition operators. A statistical ensemble can 
correctly represent the considered statistical system only when all appropriate statistical 
conditions are accurately taken into account. 

Representative statistical ensemble is a statistical ensemble that correctly represents the 
given statistical system, uniquely defining all its physical properties. The construction of 
representative ensembles for equilibrium systems is described in Ref. [27]. Then the equi- 
librium statistical operators are obtained from the conditional minimization of information 
functionals [28] . Now, we shall generalize the definition of representative ensembles to arbi- 
trary nonequilibrium statistical systems. 

Let us consider a physical system, whose correct definition requires the validity of sta- 
tistical conditions (3) for some condition operators. It is worth stressing that the latter are 
assumed to be self-adjoint, but they do not need to be compulsorily the integrals of motion. 
For instance, the normalization condition =< N > for the total number of particles 
involves the number-of-particle operator A^, which does not commute with the Hamiltonian, 
when the global gauge symmetry is broken. 

The most general way of deriving the evolution equations is by extremizing an action 
functional [29]. For a system under constraints, imposed by the statistical conditions (3), 
this implies the conditional extremization of the effective action 

A[^]= jl^m-^i^Am^ dt, (4) 

with the Lagrange multipliers Ui guaranteeing the validity of constraints (3). Defining the 
grand Hamiltonian 

H[^]=H[^]+J2^Am, (5) 
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and using Lagrangian (2), we have for the action functional (4) the form 



AVl^] ^P\x,t) i ^ ^|;{x,t) - H[^|;]^ dt . (6) 

The extremization of the effective action, defined by the variation 

5A[ijj] = (7) 
with respect to the field variables ip and ip\ implies the validity of the variational equations 



In view of action (6), this is equivalent to the equations 



(8) 



■ d , , , SHU] 

i — ipAx,t) = — r^-^ (9 

and their Hermitian conjugate. 

Thus, we obtain the evolution equations (9) for the field variables. As is seen, the 
evolution is governed by the grand Hamiltonian (5). In the Heisenberg representation, a 
field operator ■^j satisfies the Heisenberg equation of motion 

d 

i — ^j{x,t) = [ipj{x,t),H['iP]] , 

which is equivalent to Eq. (9) . Then the time evolution of the field operator is described by 
the relation 

i;J{x,t)^U+{t)^|;J{x,0)U{t) , 
with the evolution operator satisfying the Schrodinger equation 

Zj^U{t)=H[i;{x,Q)]U{t). (10) 

The evolution equations, either (9) or (10), is a necessary component for defining a 
nonequilibrium statistical ensemble. This ensemble is representative, since the evolution 
equations are derived with taking account of all statistical conditions uniquely character- 
izing the considered statistical system. The given additional constraints define the grand 
Hamiltonian (5) governing the system evolution. As is clear, the properties of a system under 
the given constraints can be essentially different from the properties of a system under other 
or without constraints. This is why the usage of representative ensembles is crucially impor- 
tant for correctly describing physical systems. The general nonequilibrium representative 
ensemble, in the case of an equilibrium system, reduces to the equilibrium representative 
ensemble with the same grand Hamiltonian (5). 
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3 Broken gauge symmetry 



Now we shall use the notion of representative ensembles, formulated above, for developing 
a correct general theory for Bose systems with broken gauge symmetry, which is associated 
with the appearance of Bose-Einstein condensate. The standard way of breaking the global 
gauge symmetry is by means of the Bogolubov shift 

^(r, t) ^(r, t) = r]{r, t) + ^i(r, t) , (11) 

in which i]{r,t) is the condensate wave function and ■?/'i(r, t) is the field operator of un- 
condensed atoms, with r being the spatial coordinate. The operators ip and -01 satisfy the 
same Bose commutation relations. The passage from to 'Xpi corresponds to a canonical 
transformation realizing nonequivalent operator representations [30]. The Fock space ^(■0), 
generated by the field operator ip^, characterizes the space of microstates for a system with- 
out Bose-Einstein condensate, while the Fock space J-'{iIji), generated by the field operator 
ipi, is the space of microstates for a system with Bose-Einstein condensate. The spaces J-'{'ip) 
and J^{ipi) are mutually orthogonal [30]. The condensate wave function ri{r,t) is the same 
as the coherent field. The field variables of condensed and uncondensed atoms are mutually 
orthogonal, such that 

lrj*{r,t)M^,t)dr^O. (12) 

The Bogolubov shift (11) is sufficient for breaking the global gauge symmetry. Note that 
the method of infinitesimal sources is not always able to break the gauge symmetry [26] . 

Having now two field variables, rj and instead of just one ip, requires to have two 
normalization conditions. One is the normalization of the condensate wave function to the 
number of condensed atoms 

No = J \v{r,t)\'dr, (13) 

where Nq, in general, can be a function of time, which is not marked explicitly just for 
brevity. Defining the operator Nq = NqI, with 1 being a unity operator in T{ipi), we can 
represent normalization (13) in the standard form of a statistical condition 

TVo = < TVo > , (14) 

as in Eq. (3). Here and in what follows, all statistical averages are accomplished over the 
space of microstates ^(V"!)- 

The second normalization condition is that one for the number of uncondensed atoms 

A^i = < ^1 > , (15) 
where the operator for the number of uncondensed atoms is 

iVi= J4{r,t)ilj,{r,t)dr. (16) 

In general, A^i can also be a function of time in a nonequilibrium system. 
The total number of atoms is 

N=No + N^, (17) 
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which could be considered as a normahzation condition, instead of Eq. (15). For two field 
variables, we have to fix not less and not more than two normalization conditions. 
The atomic densities and the related atomic fractions are denoted as 



Po = 


V ' 


Pi = 


V ' 


N 


no = 


. iVo 
N 


_ Po 
P 


rii = 


Ni pi 
N p 



where V is the system volume. For these, one has 

p = Po + Pi , no + ni 



1 . 



;i8) 



(19) 



It is assumed that po and hq are not zero in the thermodynamic limit. 

The number of condensed atoms Nq is defined so that to make the system stable. In 
equilibrium, Nq is to be found from the minimization of a thermodynamic potential. For a 
nonequilibrium system, A^o must ensure the dynamic stability of the solution for the conden- 
sate function 77(r, t) satisfying the related evolution equation. 

For a system with broken gauge symmetry, the average < ipi > may be nonzero. This, 
however, would mean that quantum numbers, such as spin or momentum, are not conserved. 
One, therefore, has to impose an additional constraint 



< Vi(r,i) > = 0. 
By introducing the condition operator 



dt , 



(20) 



(21) 



in which A(r, t) is a complex function, constraint (20) can be rewritten in the standard form 
(3) as the quantum-number conservation condition 



< A[ip] > = . 



(22) 



The grand Hamiltonian (5) has to be defined by taking into account the statistical conditions 
(14), (15), and (22), which yields 



H[r], V^i] = H[^] - poNo - PiiVi - A[V;] . 
The action functional (6), with the Bogolubov shift (11), takes the form 



(23) 



d d 
ry*(r, t) i — r){r, t) + V'I(r, t) i — V'i(r, t) 



dr - H[ri, V'l] \ dt . (24) 



The extremization of the action functional (24) implies two variational equations, for the 
condensate function, 

^%^ = 0, (25) 
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and for the field operator of uncondensed atoms, 



S4{r,t) 



(26) 



These equations, in view of the action functional (24), are equivalent to the evolution equa 
tions 

' dt^^ ' > Sr)*{r,t) 



and 



ot 5V'i(r,t) 
Let us take the energy operator in the standard form 



(27) 
(28) 



+ ^ J V^t(r, t)v;t(r', tMv - r')^(r', t)^(r, t) drdr' , 



(29) 



in which U = U{r,t) is an external field and $(r) = $(— r) is an interaction potential. To 
satisfy the conservation condition (22), the grand Hamiltonian (23) must have no linear in 
■01 terms [31]. This is achieved by chosing the Lagrange multiplier 



X{r,t) 



2m 



+ 



U + J $(r-r')|?7(r',i)|W 



(30) 



The evolution equations (27) and (28) yield 



+ /$(r-r') \7j{r',t)\'7]{r) + X{r,r') dr' 



and, respectively. 



^^^l(r,^) 



2m 



+ U -iiA V'i(r,t)+ 



(31) 



+ j $(r-r') [|77(r')|Vi(r) + r7*(r')r7(r)V'i(r') + 77(r')r?(r)V'I(r')+X(r,r')] dv' , (32) 

where, for brevity, the time-dependence is not explicitly shown in the right-hand sides of 
these equations, and the correlation operator 



X(r,r') = V^Kr')^i(r')r7(r) +^1(077(0^^1 (r) + r7*(r')V^i(r')V^i(r) + ^I(r')V^i(r')V^i(r) (33) 
is introduced. 

To obtain an equation for the condensate wave function, we need to average Eq. (31). 
For this purpose, let us define the normal density matrix 



pi(r,r') = <V'I(r')^i(r)>, 



(34) 
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the anomalous density matrix 

ai(r,r') = <^/^i(r')^i(r) > , (35) 

and the densities 

po(r) = |r/(r)p , pi(r) = pi(r,r) = < ^j(r)^i(r) > , 

(7i(r) = (7i(r,r) = <V'i(r)Vi(r)> . (36) 
The total density of atoms is the sum 

p(r)=po(r)+pi(r) (37) 

of the condensate density po(r) and of the density pi(r) of uncondensed atoms. Also, let us 
use the notation 

e(r,r') = <z^I(r')^i(r')^i(r)> . (38) 
Then the average of the correlation operator (33) becomes 

< X{r, r') > = pi(r')r7(r) + p,{v, v'^v') + a,{v, v')rj*{v') + ^(r, r') . (39) 

Finally, averaging Eq. (31), we find the equation for the condensate function 

+ I t(r-r')lp(r')»7(r)+p,(r,r')»7(r')+<7i(r,r')i)'(r')+«(r,r')] dr' , (40) 

where again, for brevity, the temporal dependence in the right-hand side is not explicitly 
shown. This is a general equation valid for an arbitrary nonequilibrium Bose-condensed 

system. 

It is worth noting that, contrary to Eq. (31), the average of Eq. (32) is not defined 
because of the following. The grand Hamiltonian does not change being complimented by 
the term 

z['ilj,] = I <C(r,i)V'I(r,i) + C(r,i)Vi(r,i)>dr, 

which is identically zero owing to constraint (20) for any complex function C(r, 0- However, 
the variation of such a term results in 



iV}(r,«) 



which, generally, is not zero. Therefore, replacing H[ri,ilJi] by H[ri, ipi] + zlipi] would lead to 
the appearance in Eq. (28) of an additional undefined term C{r,t). Equation (32) must be 
used for defining the evolution of the correlation functions < V'lV'i > < V'lV'i >■ But 
for these functions, an additive term ( does not play any role, since < CV'i >=< CV'i >— 
due to constraint (20). 
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4 Self-consistent equations 

The normal and anomalous density matrices (34) and (35), entering the evolution equation 
(40) for the condensate function, can be expressed through Green functions. One needs 
the normal and anomalous Green functions [32-34], which can be assembled in a matrix 
G(12) = [Gap{l2)] with the elements 

6-11(12) = -t< f ^i(l)^I(2) > , ^12(12) = -t< f ^i(l)^i(2) > , 

621(12) ^-i< f4{l)i;l{2) > , 622(12) = -i < f4{l)ij^{2) > , (41) 

where T is the chronological operator. Here and in what follows, the shorthand notation is 
used, denoting the set {rj,tj} just by the number j. 
We shall also need the combination 

*(123) = V'i(l)V'i(2)V'I(3) + r^(l)Vi(2)V'I(3) + V'i(l)r?(2)V'I(3)+ 

+^i(l)^i(2)r7*(3) + v{l)vi'^)4i^) + r/(l)^i(2)r/*(3) + ^i(l)r/(2)r/*(3) . (42) 

Using this, we define the binary Green function, which is a matrix 5(1234) = [i?Q,^(1234)] 
with the elements 

Sii(1234) = - < f *(123)V'I(4) > , ^12(1234) = - < f'*(123)V'i(4) > , 

^21(1234) = - < f ^+(123)^1(4) > , ^22(1234) = - < f ^+(123)^i(4) > . (43) 
We may note that 

hm < f'*(122) > = pi(2)77(l) + pi(12)?7(2) + a^{12)r)*{2) + ^(12) = < 1(12) > , 

tl2^-0 

where ti2 = ti — ^2- 

Introducing the retarded interaction potential 

$(12) = $(ri-r2)5(ii2 + 0), (44) 

we define the self-energy S(12) by the relation 

J E(13)6(32) d(3) = i J $(13)5(1332) d(3) . (45) 

Then, from Eq. (32), we obtain the matrix equation 

ri + ^ - U{1) + /xi j 6(12) - j E(13)6(32) d{3) = 5(12)i , (46) 



dti 



in which 



r = 



" 1 " 




"10" 


i = 


-1 




1 
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Equation (46) for the Green function is to be complimented by equation (40) for the con- 
densate function, which can be represented as 

' "^^^^ ^ (" S + ^^^^ ~ "^^^^ + / < T^im >] d{2) , (47) 

where Po(l) = l^(l)P- 

Equations (46) and (47) is a self-consistent set of equations derived from the extrem- 
ization of the action functional (24). Hence, these equations respect all conservation laws 
of the Hamiltonian, as is should be for the equations derived from a variational procedure 
[3,35]. It is important to stress that the self-consistency of the equations is ensured by two, 
generally different, Lagrange multipliers /^o and ^i, playing the role of the chemical poten- 
tials for condensed and uncondensed atoms, respectively. The potential Hq guarantees the 
normalization condition (14), hence, /Iq — ij,o{Nq). The normalization condition (15) defines 
A*i — f^i{Ni). The number of condensed atoms is to be such that to provide the stability of 
the system of N atoms, so that A^o = Nq{N). Because of Eq. (17), Ni = N ~ Nq. Therefore, 
A*o = A*o(-^) and Hi = Hi{N). But there is no necessity that /iq be equal to /^i, though it may 
occasionally happen. Thus, this happens in the Bogolubov approximation [21-24], which is 
valid for asymptotically weak interactions. However, for more elaborate approximations, this 
is not so. To illustrate this fact, let us consider an equilibrium system, when r]{r,t) — r]{r) 
does not depend on time. Then Eq. (47) becomes 

+ y*(r-r') Wr')i)(r)+pi(r,r')i;(r')+<Ti(r,r')i)*(r')+«r,r')l dr' = /i„i,(r) . (48) 
The grand thermodynamic potential is 

Q= -rinTrexp(-/5i/[r;,^i]) , (49) 

and the free energy being 

F = -TlnTrexp{-/3(iy[t/'] - A[t/'])} , (50) 

where P — T'^ is inverse temperature. For potentials (49) and (50), we have the relation 

n = F- HoNo - /xiiVi . (51) 

The number of condensed atoms is defined by the stability condition implying the minimum 
of the grand potential (49), 

More generally, the condensate wave function, assuring the system stability, is defined by 
the condition 

6^) = ^^^''^^- ^^^^ 
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The latter, owing to Eq. (27), gives Eq. (48). 

Since A^o = Nq{N) and Ni = Ni{N), it is the total number of atoms N that can only be 
fixed for a system, which requires the relation 

n^F-iiN , (54) 

in which /j, is the system chemical potential. Comparing Eqs. (51) and (54), immediately 
results in the definition of the system chemical potential 

II = no/iQ + ni/ii . (55) 

For the latter, one has the standard thermodynamic equations 

Let us emphasize that, since Nq and Ni are uniquely defined through the total number 
of atoms A^, neither Nq nor Ni can be treated as independent thermodynamic variables, 
because of which, generally, ^ f^i- 

To give a more explicit demonstration that hq ^ /xi, let us turn to a uniform system, 
when the density matrices pi(r, r'), (Ti(r, r'), as well as ^(r, r'), depend solely on the difference 
r — r'. Then the densities p\ = pi(r, r). o"i(r. r), and ^i(r, r) are constants, together with the 
order parameter rj — r]{r). Also, normalization (14) yields rj — ^/po. From Eq. (48), we find 

A(o = p$o+ /(nik + (7fe + 4^)$fc , (56) 

where n^, a^, and are the Fourier transforms of the corresponding quantities $(r), 
Pi(r, 0), cri(r, 0), and ^(r, 0), respectively. For the Fourier transforms ^^/^(k, uj) of the Green 
functions Gap{12), we find the former from Eq. (46). The poles of Ga/3(k, a;) define the 
single-atom spectrum, which is gapless, provided that 

/xi = En(0,0)-Si2(0,0) , (57) 

where Ea/3(k, a;) is the Fourier transform of the self-energy Eq,/3(12). Note that the gapless 
spectrum is a necessary requirement for the existence of a stable system with Bose-Einstein 
condensate [30]. Equation (57) is the Hugenholtz-Pines relation that could be derived either 
from thermodynamic equations [33] or from the Ward identities with respect to the variation 
of gauge [24] . 

One often considers Bose-condensed systems in the contact-potential approximation, 
when 

$(r) = $o<^(r) f $0 = 47r ^) , 
with as being a scattering length. Then — and Eq. (56) gives 

//o= (p + Pi + (7i + ^) $0 , (58) 
V VPoJ 
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where 

f dk f dk f dk 

^^=y^'=(2^' ^^=y^'^(2^' ^^y^^(2^- 

Clearly, there is no such a general law that would require the identity of /iq and /^i. 

To be even more specific, let us consider the Hartree-Fock-Bogolubov approximation, 
when ^(r, r') = 0. Then for Eq. (57), we find 

f dk 
/xi = p$o + J {rik - crk)^k ■ (59) 

In the case of the contact potential, we obtain 

^0 = (P + Pl + (7i)$o , //i = (p + pi - (7i)$o ■ (60) 

As is evident, /xq 7^ A^i- 

In the standard approach, without using the representative ensemble, one assumes that 
/Iq equals ^i, hence, according to Eq. (55), equals But then one comes to the explicit 
inconsistency in Eqs. (60). A common way of trying to treat this inconsistency is by 
neglecting the anomalous average ai, calling this the "Popov approximation". However, this 
unjustified trick has nothing to do with Popov, as is easy to infer from his original works 
[36-38]. Moreover, it is easy to show that at low temperatures the anomalous average ai can 
be much larger than the normal average pi, that is, neglecting ai has nothing to do with a 
reasonable approximation [30,39]. 

By employing the representative ensemble, as is done in the present paper, we introduce 
two Lagrange multipliers, and fii, for two normalization conditions (14) and (15). These 
multipliers are not obliged to be equal, but are to be such that to render the whole theory 
completely self-consistent. 



5 Topological coherent modes 

Equation (48), defining the condensate function 77(r), corresponds to a stable equilibrium 
system, when condensing atoms pile down to the ground-state quantum level. This equation 
can be generalized for describing arbitrary stationary states ?7n(r), enumerated by a quantum 
multi-index n, so that, for the contact potential, we have 



2m 



+ U{r) 



+$o{[|^„(r)r + 2pi")(r)]r7„(r) + ai")(r)r7:(r) + e("nr)} = Kr7n(r) , (61) 

where p^^\r), ai^\r), and ^'•"•'(r) arc the solutions to the self-consistent system of equations 
for pi(r), cri(r), and ^(r) = ^(r, r), respectively, under the condition that ?7(r) is replaced by 
r]n{r). The ground-state energy level corresponds to 



Eq = min ^ po , 



(62) 
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with Tjolr) = r){r). But, in general, there can exist a whole set of coherent states rin{r), with 
a spectrum of energies E^- The condensate functions rjnir) will be called the topological 

coherent m,odes. 

It is clear that, for a system at absolute equilibrium, the sole pertinent mode is the 
ground state 77(r) . In order to generate excited coherent modes, it is necessary to deal with 
a nonequilibrium system. Then one has to consider the time-dependent equation (40) for 
the condensate function, which, for the contact potential, reads as 

d / \ 

+$0 {[po(r, t) + 2pi(r, t)] ^(r, t) + ai(r, t)'n\v, t) +i{v,t)} , (63) 

where U = U{r, t) is a time-dependent external potential. The latter can be represented as 
a sum 

C/(r,t) = f/(r) + \/(r,t) (64) 

of a trapping potential f/(r) and a modulating potential V{r,t). 

If in Eq. (63) we set pi(r,t) = 0, cri(r,t) = 0, and $,{r,t) = 0, then we come to the 
usual temporal Gross-Pitaevskii equation [1] . The generation of various topological coherent 
modes on the basis of the latter equation has been investigated earlier [6-20,40]. However, 
in general, we have to analyse the full Eqs. (61) and (63). 

Suppose that at the initial time t = 0, the atomic system is in equilibrium, so that 

77(r,0) = 77o(r) = 77(r). (65) 

Assume that we wish to generate a nonground-state condensate, corresponding to a topo- 
logical coherent mode labelled by the index n — rii, with the energy Ei = En-^. For this 
purpose, we impose an alternating modulating potential 

V(r, t) = Vi(r) cosa;t + V2{r) sin cut , (66) 

with the frequency cu tuned close to the transition frequency cui = Ei — Eq, which implies 
the resonance condition 

— < 1 {Au = 00-001) . (67) 
00 

Imposing an external alternating field will, of course, destroy the ground-state condensate 
by means of two processes. One is the resonant process of transferring condensed atoms from 
the ground state mode, with the energy Eq = ^q, to the chosen excited mode, with the energy 
El. And also, there will be nonresonant processes of taking atoms to other nonresonant 
coherent modes, as well as the process of transferring atoms from the condensate to the 
cloud of noncondensed atoms. Estimates show [15] that it is feasible to arrange such a setup 
that the resonant generation would occur much faster than other nonresonant processes. 
Assuming this, we may consider the solution of Eq. (63) at times shorter than the critical 
time tc, when the resonant generation prevails, and the number of condensed atoms stays 
practically constant. 



J \r]{r,t)\'^dr = J \r]{r)\'^dr = No . {6i 
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The topological coherent modes are also normahzed to A^o, because of which we can set 

rin{r) = 0Vo (^„(r) , J \ipn{r)\^dr = 1 . (69) 
Then we look for the solution of Eq. (63) in the form of the mode expansion 

7j{r,t) = J2Cn{t)vn{r)e-'^"\ (70) 

n 

in which ujn = En — Eq and the coefficient functions C„(t) are assumed to be slow functions 
of time, as compared to the exponential exp(— ia;„t), which means the condition 



1 



dCri 



dt 



< 1 . 



(71) 



From the normalization condition (68), with expansion (70), after averaging over time, while 
keeping the quasi- invariants Cn{t) fixed, we get 



(72) 



Substituting expansion (70) into Eq. (63), we multiply the latter by 77*(r) exp(ia;„i), 
integrate over r, and average over time, following the averaging technique [6,14,15]. To this 
end, we need the following integrals 



^oNo j Wm{r)? [2|<^n(r)P - km(r)|^] dv , 



Pmn = j ¥^;(r) mv) - miv)] (^„(r) dv , 



olnnit) ^ - $0 / ^^^(r) | 2 [pt\v) - pi(r, t)] <pn{v) + <jt\v)<pl{v) + 2^ I dv , 



^qNq j \ipn{v)\^dv , = Aa; + aoo - ctn . 



We presume that, at the resonant stage, when t <^ tc, the quantity Q;„„(t) is a slow function 
of time, such that 



tr 



dar 



dt 



< 1 . 



Then, introducing the functions 

Cn{t) = Cn{t) exp {iann{t)t} , 

we, finally, arrive at the system of two equations 

• dCo I 1 2 ^ n I ■ \ \ 

% — = q;oi|ci| Co + - /^oiCi exp(^Aoit) , 



(73) 



(74) 



dt 

Icq 
dt 



■ I |2 ( ■ \ ^\ 

I — = aio|co| ci + - iJoiCocxp{-iAoit) . 



(75) 
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Equations (75) have the same form as studied in the previous pubhcations [6,8,10,14-16,20] 
describing the resonant generation of topological coherent modes on the basis of the Gross- 
Pitaevskii equation. 

The reduction of the general equation for the condensate function (63) to the system of 
Eqs. (75) demonstrates that the topological coherent modes, corresponding to nonground- 
state condensates, can be resonantly generated also in the presence of the cloud of uncon- 
densed atoms. The procedure of the resonant generation requires the validity of several 
conditions discussed above. In particular, the resonant process must be much faster than 
nonresonant ones. After the critical time tc, power broadening destroys the resonant proce- 
dure and nonresonant processes become prevailing. According to estimates [15], the critical 
time is of the order of tc ~ Lojip? + /3^), where a = (aoi + cno)/2 and P = \Poi\. At the 
initial stage of time t <^ t^, the resonant generation is feasible. The critical time tc can be 
made comparable to the lifetime of atoms in a trap [15]. 
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